Understanding cuprate superconductors with spontaneous nodal gap generation 
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We study the spontaneous gap generation for gapless nodal fermions within an effective gauge 
field theory of high temperature superconductors. When superconductivity appears, the gauge 
boson acquires a finite mass via Anderson-Higgs mechanism. Spontaneous nodal gap generation 
takes place if the gauge boson mass £ is zero or less than a critical value £ c but is suppressed 
by a larger gauge boson mass. The generated nodal gap prohibits the appearance of low-energy 
fermion excitations and leads to antiferromagnetic order. Using the fact that gauge boson mass £ is 
proportional to superfluid density and doping concentration, we build one mechanism that provides 
a unified understanding of the finite single particle gap along the nodal direction in lightly doped 
cuprates, the competition and coexistence of antiferromagnetism and superconductivity, and the 
thermal metal-to-insulator transition from the superconducting state to the field-induced normal 
state in underdoped cuprates. 
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I. INTRODUCTION 



Recently, angle-resolved photoemission spectroscopy 
(ARPES) measurements were performed on lightly doped 
cuprates and found a full gap over the whole Brillouin 
zonei. This finding is unexpected since the d-wave sym- 
metry of the gap/psedogap of high temperature super- 
conductors has widely accepted^ and previous extensive 
ARPES measurements always revealed gapless excita- 
tions in the (±7r/2, ±7r/2) direction^. The nodal gap 
decreases upon doping and disappears when supercon- 
ductivity emerges as the ground state. In some materials 
with superconductivity being its ground state, the nodal 
gap is also observed. This indicates that there seems to 
be a competition and possible coexistence between the 
nodal gap and superconductivity. This behavior is quite 
similar to the evolution of orders with doping concentra- 
tion. Antiferromagnetism is the ground state of undoped 
and lightly doped cuprate superconductors, but disap- 
pears when superconductivity emerges, thus there is a 
competition between antiferromagnetism and supercon- 
ductivity. Because of this competition, when supercon- 
ductivity is suppressed, say by magnetic fields, antifer- 
romagnetism has a chance to appear in a superconduc- 
tor. The field-induced local antiferromagnetic order has 
been confirmed by recent neutron scattering and scan- 
ning tunnel microscopic (STM) experiments^*^. In 
particular, antiferromagnetic order is found not only in 
the Abrikosov vortices but also in the superconducting 
region around the vortices. This finding strongly sup- 
ports the coexistence of antiferromagnetism and super- 
conductivity. On the other hand, extensive heat trans- 
port measurements show that a residual linear term ex- 
ists in the whole superconducting region at T — > 0, but it 
decreases with decreasing doping concentration and ap- 
proaches zero as superconductivity disappears&2*i2*ii*i£. 
When an external magnetic field is present perpendicular 



to the Cu02 plane, the thermal conductivity of under- 
doped cuprates decreases with increasing magnetic field 
and finally vanishes when the magnetic field is beyond 
the up critical field H. c %. Thus cuprate superconductors 
exhibit a thermal metal-to-insulator transition upon go- 
ing from the superconducting state to the field-induced 
normal stated. 

We believe that the above phenomena are universal to 
all high temperature supercodncutors and more impor- 
tantly they are governed by the same physical mecha- 
nism. In this paper, we argue that all these experimental 
results can be understood if the gapless nodal fermions 
acquire a finite gap. Such a spontaneous gap generation 
for nodal fermions is achieved by coupling the fermions 
to a gauge field which naturally appears as a result of 
strong electron correlation effect when we go beyond the 
slave-boson mean-field treatment of t-J model. 

Spin-charge separation and emergent gauge fluctua- 
tion are two key concepts in our scenario. When spin 
and charge degrees of freedom are separated, the ex- 
citations are spin-carrying spinons and charge-carrying 
holons rather than ordinary electrons. The pairing of 
spinons is responsible for the observed <i-wave energy 
gap/pseudogap. The d-wave gap vanishes at the nodes, 
so the low-energy fermion excitations are effectively gap- 
less and hence can be described by relativistic massless 
Dirac fermions. Superconductivity is realized once the 
holons undergo Bose condensation at low temperatures. 
The spinons are connected to the holons due to the ex- 
change of an emergent gauge field although they do not 
interact directly. Thus, the low-energy behavior is domi- 
nated by an interacting system consisting of gapless Dirac 
fermions, holons and an emergent gauge field. 

A finite nodal gap is generated spontaneously when 
the gauge field binds gapless nodal fermions into stable 
fermion- anti- fermion pairs. Exciting single fermions from 
the nodal direction of the Brillouin zone requires a finite 
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energy cost which is responsible for the finite nodal gap 
observed in ARPES experiments. Since the single par- 
ticle spectrum is fully gapped in the whole momentum 
region, no free fermions can exist at low temperatures 
and consequently there should not be a linear term for 
the thermal conductivity, which is present in a pure d- 
wave superconductor with gap nodes. Spontaneous nodal 
gap generation can take place if the holons are absent 
or free. When superconductivity emerges as the ground 
state, the local gauge symmetry is broken by the holon 
condensation and the gauge boson becomes massive via 
Anderson-Higgs mechanism. We found a critical value 
for the gauge boson mass. Spontaneous nodal gap gen- 
eration takes place if the gauge boson mass is less than 
the critical value but is suppressed when the gauge boson 
mass becomes larger than the critical value. Then there is 
a competition between spontaneous nodal gap generation 
and superconductivity. The linear term for thermal con- 
ductivity in the superconducting state observed by heat 
transport measurements reflects the suppression of spon- 
taneous nodal gap generation by superconductivity. An- 
other important consequence of spontaneous nodal gap 
generation is that it enhances the antiferromagnetic spin 
correlation greatly and actually leads to long-range or- 
der. Thus the antiferromagnetic order should have the 
same doping dependence with the nodal fermion gap: it is 
present in low doping region and disappears in high dop- 
ing region after superconductivity emerges. This is ex- 
actly what happens in cuprates superconductors. Since 
the gauge boson mass is proportional to the superfluid 
density, the coexistence of spontaneous nodal gap gen- 
eration and a small gauge boson mass leads to an very 
important result that the antiferromagnetic order can co- 
exist with superconductivity when the superfluid density 
is less than the critical value. Thus the spontaneous 
nodal gap generation provides a unified explanation of 
the experimental results mentioned at the beginning of 
this paper. 

This paper is organized as follows: In Sec. (II), we first 
give a brief review on spin-charge separation, slavo-boson 
treatment of t-J model and the gauge theory approach 
to cuprate superconductors, then we discuss chiral sym- 
metry breaking in the presence of massive gauge boson 
and estimate the critical gauge boson mass which sepa- 
rates chiral symmetric and symmetry broken phases. In 
Sec. (Ill), we discuss the explanation of the ARPES ex- 
periments, the heat transport behavior and the relation- 
ship between antiferromagnetism and superconductivity 
respectively. The paper ends with a summary. 



II. DYNAMICAL GAP GENERATION FOR 
GAPLESS NODAL FERMIONS 

Shortly after the discovery of high temperature super- 
conductors, it was correctly recognized that these ma- 
terials are doped Mott insulators and hence can not be 
properly understood without considering the strong cor- 



relation effect. Anderson^ 3 , proposed that due to quan- 
tum fluctuations the ground state of undoped cuprates is 
more probably some kind of quantum liquid of spin sin- 
glets, called resonating valence bond (RVB) state. The 
investigation of RVB state was carried out within the 
t — J model which is derived from the more general three- 
band Hubbard modeUi The strong correlation nature of 
cuprate superconductors is reflected in a no-double occu- 
pancy constraint which says that there is no more than 
one electron at one lattice site due to the strong Coulomb 
repulsion between electrons. By decomposing the elec- 
tron operator c\ a to the product of a spinon operator 
j\ a (neutral fermion) and a holon operator hi (charged 
boson) 

cl = flh, (1) 
the no-double occupancy constraint can be written as 

£/L/w + &l&i = i (2) 



which is easier to be treated analytically. This decompo- 
sition is the crucial step underlying the so-called slave- 
boson mean field treatment of t — J model. A four- 
fermion interaction term appears in the t-J model af- 
ter replacing electron operators with spinon operators 
and holon operators. This term can be treated by in- 
troducing three order parameters Xij — (fia-ficr)i ^ij = 

(filfjl - filfjt)> and Vij = b\bj. A phase diagram can 
be obtained based on this mean field approachi^i^. The 
flux phase was found to be locally stable^ and very in- 
teresting for its applicability to cuprate superconductors. 
Due to its d-w&ve spinon gap symmetry, the low-energy 
excitations are actually gapless Dirac fermionsiii^. The 
quantum fluctuations around this mean field state in- 
cludes a massless U(l) gauge field. The low-energy ef- 
fective behavior of cuprate superconductors thus can be 
well described by the three-dimensional quantum electro- 
dynamics (QED3). 

Although the U(l) formulation captures some impor- 
tant properties of high temperature superconductors, it 
was shown by Wen and Leei^ not to connect smoothly 
to the hall-filling material in which an exact SU(2) lo- 
cal symmetry was found. An SU(2) treatment of the 2D 
t-J model was constructed to describe the physics of un- 
doped and underdoped cuprates in a unified wayiS. This 
treatment gives rise to several mean-field phase diagrams 
among which the staggered flux phase is applicable to 
cuprate superconductors. Taking quantum fluctuations 
into account leads to a low-energy effective theory 2 ^ that 
consists of a massless U(l) gauge field, a two-component 
bosons, and a massless fermions excited from the gap 
nodes of d-wave spinon pairs. 
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A. Dynamical chiral symmetry breaking 

If the fermions are massless the theory respects a chiral 
symmetry. However, the QED 3 theory has a rather pe- 
culiar property that the massless fermions can acquire a 
dynamically generated mass when its flavor is less than a 
critical value. The mass term of fermions spontaneously 
breaks the chiral symmetry which leads to a massless 
Goldstone boson according to the Goldstone theorem. 

At half-filling, the low energy physics is dominated by 
the interaction of nodal fermions with U(l) gauge fields^ 



N 



(dp - ia M )7 M i/v- 



(3) 



The Fermi field is a 4 x 1 spinor representing the 
gapless nodal fermions. The 4 x 4 7 M matrices obey the 
algebra, {7^, 7„} = 25^, and for simplicity, we let v CT . p = 
1 (fj,, v = 0, 1, 2). At half-filling case, there are no holons 
or the holons are all confined because of the presence of 
a very large charge gap. 

Chiral symmetry breaking is a nonperturbative phe- 
nomenon and can not be obtained within any finite or- 
der of the perturbation expansion. The standard ap- 
proach to this problem is to solve the self-consistent 
Dyson-Schwinger (DS) equation for the fermion self- 
energy. The inverse fermion propagator is written as 
S~ 1 (p) = ij-pA (p 2 ) +X (p 2 ) , A(p 2 ) is the wave-function 
renormalization and S(p 2 ) the fermion self-energy. The 
propagator of a massless fermion is simply S^ 1 (p) = vy-p. 
A(p 2 ) and £(p 2 ) appear due to the renormalization effect 
caused by interaction with gauge field. The self-energy 
function £(p 2 ) represents the interaction induced fermion 
mass and is determined by a set of DS integral equations. 
If the DS equation for S(p 2 ) has only vanishing solu- 
tions, the fermions remain gapless and the Lagrangian 
respects the chiral symmetries ip — > exp(i0j3^)ip, with 
73 and 75 two 4x4 matrices that anticommute with 7^ 
(p = 0,1,2). If the DS equation for T,(p 2 ) develops a 
squarely integrable nontrivial solution^!, then the orig- 
inally massless fermions acquire a finite mass. The DS 
equation is extremely complicated and hence can never 
be treated without making proper approximations. To 
the lowest ordes 2 ^ we take A = 1 neglecting all higher- 
order corrections and approximate the vertex function by 
the bare 7^. After making these approximations, the DS 
equation has the form 



d 3 k yi*D liU (p-k)'£(k 2 )'f 



(2tt) 3 



k 2 + E 2 (fc 2 ) 



(4) 



Appelquist et al. showe d 22,2 ? that this DS equation has 
solutions only for N < N c = 32 /it 2 . Since the physical 
fermion flavor is 2 representing the number of spin com- 
ponent, the massless fermions actually acquire a finite 
mass. This mechanism is called dynamical chiral symme- 
try breaking and has been studied for many years in par- 
ticle physics as a possible mechanism to generate fermion 



mass without introducing annoying Higgs bosons. The 
presence of a critical fermion flavor can be understood as 
follows. The effective coupling of gauge field is propor- 
tional to 1/N. It is strong enough to form fermion pairs 
when the physical fermion flavor is less than N c , while 
for large N the coupling is too weak. 



B. Effect of gauge boson mass on chiral structure 

When holes are doped into the CuO plane, the holons 
are excited and hence an additional coupling between 
holons and the gauge field appears. We have shown24 
in a gauge invariant way that the gapless fermions can 
acquire a finite gap when the additonal holons are not 
Bose condensed. When superconductivity emerges, the 
gauge boson becomes massive. In order to understand 
the properties of nodal fermions we need to investigate 
the effect of the gauge boson mass on chiral symmetry 
breaking. 

We write the action of bosons in the standard 
Ginsberg-Landau form 



4m 6 



\{d,-ia,)b\ 2 - a \b\ 2 -l-\b\\ (5) 



Note that this is not the popular model that has been 
extensively studied in the literature 2 ^. In previous treat- 
ment, a non-relativistic model has been used to describe 
the interaction of holons and gauge fields. In such a 
model, the density fluctuations of holons screen the tem- 
poral component of the gauge field which becomes mas- 
sive and hence is ignored. However, this treatment de- 
stroys the gauge invariance of the theory: the result ob- 
tained in the Landau gauge is qualitatively different from 
that in the Feynman gauge 2 ^. This inconsistency might 
be a result of using an inappropriate Lagrangian for the 
holons. At present, it is not possible to derive an effec- 
tive action for the holons rigorously. In this paper, we 
use the relativistic scalar QED to describe the holons be- 
cause it is the most general field theoretic model for a 
scalar field and it can lead to a gauge invariant critical 
fermion flavor. 

In the Lagrangian, (3 is always positive while a can be 
positive or negative corresponding to normal and super- 
conducting states, respectively. For a > 0, the holons are 
free and the Lagrangian Lb is invariant under the local 
U(l) gauge transformation 



(6) 



(7) 



where 9 (x) is an arbitrary function. When the holons 
undergo Bose condensation, a becomes negative and the 
ground state occurs at 



(b) 



(8) 
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Thus, the local gauge symmetry is spontaneous broken 
due to ground state degeneracy. We could write the holon 
field as 

b(x)=b e ie W, (9) 

where 9 is the phase of the order parameter which is just 
the gapless Goldstone mode associated with the sponta- 
neous gauge symmetry breaking. The appearance of this 
gapless mode used to a big puzzle to physicists since no 
such modes had been observed in superconductors. This 
inconsistency could be eliminated by the Anderson-Higgs 
mechanism-^. Its essential idea is to perform the follow- 
ing gauge transformation 

a M — > a M + dfj.9, (10) 



Let q be the gauge boson momentum, then q 2 = (p — 
k) 2 = p 2 + k 2 — 2pkcos9. Tt(q 2 ) is the vacuum polariza- 
tion of the gauge boson, which is originally introduced 
to overcome the infrared divergence. In the present case, 
the gauge field has no kinetic energy term and its dy- 
namics comes from integrating out matter fields. As a 
result, only ir(q 2 ) appear in Dr(q 2 ), so 

D T 1 (q 2 )=q 2 7T(q 2 )+e- (16) 

The vacuum polarization consists of two parts corre- 
sponding to fermion contribution 7Tpand holon contribu- 
tion itb respectively, which are 

W) = ^, (17) 



taking advantage of gauge freedom of the theory. Then 
the Goldstone mode 9 is removed and there appears a 
term as follows 



— b 2 a^. 



(11) 



It is easy to see that the gauge boson now acquires a 
finite mass 



£ = — Ps 
mb 



(12) 



after absorbing the gapless Goldstone mode. This is the 
famous Anderson-Higgs mechanism which when general- 
ized to non-Abelian gauge theories constitutes the foun- 
dation of the electro-weak Standard Model. The physical 
meaning of the finite gauge boson mass can be seen from 
the equation for magnetic field 



V 2 B = £B. 



(13) 



Comparing this with the standard London equation, we 
know that £ = A^ 2 with A^ the London penetration 
depth. The solution of this equation is an exponen- 
tially damping function indicating that the magnetic field 
can penetrate the superconductors only to a finite depth 
equal to the inverse gauge boson mass. 

After the gauge boson acquires a finite mass, its cou- 
pling strength is weakened and might not be able to form 
fermion pairs. To testify the correctness of this naive ex- 
pectation, we now investigate the effect of a finite gauge 
boson mass on chiral symmetry breaking by studying 
the DS equation with massive gauge boson propagator 
and seek the critical coupling constant. The gauge boson 
propagator in Landau gauge is 



ZV(q) = D T (q 2 ) 



with 



D T \q 2 ) = q 2 [l + A<l 2 )]+e 



(14) 



(15) 



1 



(18) 



to the one-loop approximation (Note that the holon mass 
is ignored in the polarization ttb since it does not affect 
the chiral structure 2 ^). Adding them up leads to the total 
vacuum polarization 



N + 1 



Then we have 



D^(q 2 ) = qMq 2 )+e 



N+l 



(q + v), 



with 



V 



N + l 



(19) 



(20) 



(21) 



Now we can write down the propagator of gauge boson 
explicitly 



(N + l)(\q\+r)) 



5uu 



q 2 



(22) 



We now substitute this expression into Eq. (4). After 
performing angular integration and introducing an ultra- 
violet cutoff A we have 



S(p 2 ) = A / 
Jo 



dk 



fc£(fc 2 ) 



x- [p 

P 



k 2 
-k 



Y?{k 2 ) 
\p-k\- 



77 In 



\p-k\+V 



3) 



where A = 4/ (A + l)7r 2 serves as an effective coupling 
constant. 

This integral equation can be investigated by bifurca- 
tion theory and parameter imbedding method. In order 
to obtain the bifurcation points we need only to find the 
eigenvalues of the associated Frechet derivative of the 
nonlinear DS equation^. Those eigenvalues that have 



5 



odd multiplicity are the bifurcation points. The first bi- 
furcation point is just the critical coupling strength at 
which a nontrivial solution of the DS equation develops 26 . 
Once we obtained the critical coupling constant, then we 
can get the critical fermion flavor that separates the chi- 
ral symmetry breaking phase and the chiral symmetric 
phase. 

Numerical calculations^ 6 , found that the critical 
fermion number is a monotonically increasing function 
of A/77. It conforms the naive expectation that a fi- 
nite mass of the gauge boson is repulsive to gap gen- 
eration for fermions. For small A/77 the critical number 
is smaller than physical fermion number 2, so fermions 
remain gapless. When A/77 increases, the critical number 
increases accordingly and finally approaches a constant 
value larger than 2. Thus we can conclude that the spon- 
taneous gap generation for fermions takes place when the 
gauge boson mass is zero and very small but is destroyed 
when the gauge boson mass is larger than a critical value. 
Including the wave function renormalization A(p 2 ) does 
not change this behavior, but changes the critical value 
quantitatively. We have performed calculations in an ap- 
propriate non-local gauge after taking A{p 2 ) into account 
and found that ^ is between 10~ A eV and 10~ 3 eK Al- 
though it is hard to determine its exact value, we will 
show that this uncertainty only has minor influence on 
the critical values of observable quantities such as the 
superfluid density and the doping concentration. 

We know from these results that there is a competi- 
tion between two kinds of spontaneous symmetry break- 
ing: chiral symmetry breaking and spontaneous gauge 
symmetry breaking. Chiral symmetry breaking is associ- 
ated with the mass generation of massless fermions, while 
spontaneous gauge symmetry breaking is caused by holon 
condensation and generates a mass for the gauge boson. 
If the gauge boson mass is less than the critical value 
but nonzero, there is a coexistence of chiral symmetry 
breaking and spontaneous gauge symmetry breaking. 



C. Critical point between chiral symmetric and 
symmetry breaking phases 

We now would like to discuss the critical value for the 
gauge boson mass since it plays a crucial role in deter- 
mining the transition between chiral symmetric and sym- 
metry breaking phases. The gauge boson mass is not a 
physical quantity that can be observed by experiments. 
In order to explain experiments with our mechanism, we 
should make a connection between the gauge boson mass 
with the superfluid density and the doping concentration, 
which arc important physical quantities in describing the 
superconductors. 

The relationship between gauge boson mass with su- 
perfluid density we obtained in Eq. (12) is very useful in 
describing the effect of superconducting condensation on 
spontaneous gap generation for the nodal quasiparticles. 
Now we wish to relate the gauge boson mass to dop- 



ing concentration S. Careful London penetration depth 
measurements and optical conductivity experiments 2 ^ 
showed that the superfluid density is proportional to the 
doping concentration, that is 

Ps (T = 0) = ^, (24) 

where a is the lattice constant. The optical conductivity 
experiments reflects the response of the system to electro- 
magnetic field, instead of the internal gauge field. This 
seems to be an obstacle to use the above equation. How- 
ever, we expect it works well for both the internal gauge 
field and the electromagnetic field since they have the 
same gauge structure. Thus, we obtain the relationship 
between gauge boson mass and doping concentration 

5 = m b a 2 £. (25) 

Based on this formula, the doping dependence of many 
physical quantities can be described by their dependence 
of the gauge boson mass. Note that this property is spe- 
cial to cuprate superconductors which are believed to be 
doped Mott insulators. 

Next we would like to calculate the critical value of the 
doping concentration which separates the chiral symmet- 
ric and symmetry breaking phases. As will be seen in the 
following discussions, this is an important quantity in un- 
derstanding experiments. Our calculation of DS equation 
has given the critical gauge boson mass. From the nu- 
merical results in the nonlocal gauged, we know that ^ 
is between 10 _4 eV and 10~ 3 ey. The continuum U(l) 
gauge theory is the effective low-energy theory of high 
temperature superconductors, hence the lattice provides 
a natural ultraviolet cutoff. However, we can also choose 
a = (N + l)/8 as the ultraviolet cutoff since in QED3 
all physical quantities damps rapidly above this energy 
scaled 3 -. Remember that we have defined 77 = 8£ 2 / (JV+1). 
Putting all these together, we get the critical doping con- 
centration 



The mass of the holons is determined by the hopping in- 
tegral t in the t-J model. In the tight-binding treatment, 
we have 

(2m b a 2 y 1 = t h = 0.122ey, (27) 

which was obtained by Lee and Wen 28 in the case of 
YBC0 6 . 95 . If ^ = 10" 4 eV then the critical doping 
concentration is 8 C = 0.03; while if ^ = 10~ 3 eV, then 
5 C = 0.05. Since the antiferromagnetic order disappears 
generally at 0.03, the critical 5 C we obtained is in good 
agreement with experiments and the inability in deter- 
mining the exact value of ^ does not affect the reliability 
of our conclusion. 

Spontaneous gap generation takes place for doping 
concentration less than <5 C , no matter the ground state is 
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superconducting or not. For doping concentration larger 
than S c , if the ground state is not superconductivity, then 
spontaneous gap generation also takes place; however, if 
the ground state is superconducting, then spontaneous 
nodal gap generation is suppressed by superconductiv- 
ity. In the superconducting state with doping concen- 
tration higher than 8 C , there is possibility that the su- 
perfluid density is reduced to below its critical value p sc 
by some means other than decreasing doping. If this re- 
ally happens, then spontaneous gap generation can also 
take place. This coexistence of superconductivity and 
spontaneous gap generation is very important and will 
be discusses in the next section. 



III. EXPLAINING EXPERIMENTS WITH 
SPONTANEOUS NODAL GAP GENERATION 

Since we know the effect of holons and gauge boson 
mass on the gap structure of nodal fcrmions, we are 
equipped for understanding the exotic experimental find- 
ings we have mentioned in the Introduction. These ex- 
periments include the finite single particle gap along the 
nodal direction, the low temperature quasiparticle heat 
transport, the doping dependence of antiferromagnetism 
and its relationship with superconductivity. The com- 
mon feature of these phenomena is that they are all dom- 
inated by the behavior of nodal quasiparticles. It is thus 
not surprised that they can be understood by one sin- 
gle physical mechanism. Although our mechanism is still 
rather qualitative, its efficiency in explaining these exper- 
imental results gives us confidence that it does capture 
some essential physics of hign temperature superconduc- 
tors. 



A. Finite energy gap in nodal direction 

Now we discuss the application of our theory to 
ARPES experiments. In RVB theories, all the single par- 
ticle gaps are caused by the pair formation of spinons. 
The spinon gap exists at any doping concentration that 
is less than optimal doping and has d-wave symmetry. 
ARPES measurements play an significant role in study- 
ing the gap symmetry since it is momentum dependent. 
Most previous ARPES measurements have been limited 
to optimally doped Bi2Sr2CaCu208+a; because these ma- 
terials have ideal surfaces which are required by ARPES 
measurements. However, to understand high tempera- 
ture superconductivity, it is necessary to know the elec- 
tronic structure of underdoped, lightly doped and un- 
doped cuprates. 

Recently, ARPES measurements have been per- 
formed in the lightly doped cuprates including hole- 
doped La2_ a; Sr 2: Cu04 and Ca2- K Na a; Cu02Cl2 and 
electron-doped Nd 2 _ a; Ce :I ;Cu04i. The spectra from 
Ca2_2,Na x Cu02Cl2 reveals that a finite gap exists in the 
nodal direction at x = 0.05 and becomes smaller with 



increasing doping concentration. It closes at x = 0.12 
corresponding to a critical temperature T c = 22K. For 
La2- a; Sr a ;Cu04, a clear nodal gap is observed at doping 
concentrations x = 0.01 and x = 0.02, but it closes at x = 
0.03 which is less than the critical doping x c = 0.05 where 
superconductivity starts to emerge as the ground state. 
In the case of slightly electron-doped Nd2-xCea;Cu04, 
the gap in the nodes is observed at doping x — 0.025 and 
x = 0.04, and it does not exist at doping x = 0.08 and 
x = 0.10, well below the critical doping concentration 
x c = 0.12. Such a doping dependence of the excitation 
gap in the nodal direction is observed in three different 
cuprate superconductors and hence should be a universal 
phenomena. These observations are rather striking be- 
cause based on the d-wave symmetry of gap/pseudogap, 
the energy gap should vanish along the nodal directions. 

Although the nodal gap in the three cuprates has the 
same dependence of the doping concentration, the in- 
consistency of critical doping at which the nodal energy 
gap closes and critical doping at which superconductivity 
emerges brings a difficulty to understand the relationship 
between the nodal gap and high temperature supercon- 
ductivity. However, the difficulty is not that severe as it 
seems to be, because all the above ARPES measurements 
are performed at a finite temperature T — 15K rather 
than nearly zero temperature. It is very probable that 
the true ground state of all the lightly doped cuprates 
has the same electronic structure and the above ARPES 
observed inconsistency is caused by thermal fluctuations 
which arc different from material to material. 

We speculate that this picture is what really happens 
in high temperature superconductors. In the spirit of 
our results, at doping concentration that is lower than 
the critical doping at which superconductivity appears, 
the instability caused by strong gauge fluctuations always 
generates a finite gap for the nodal fermions. When su- 
perconductivity emerges, if its superfluid density is less 
than some critical value, the nodal fermions also have a 
finite gap. This finite gap is suppressed by a larger dop- 
ing concentration, i.e., a larger superfluid density. This is 
a universal picture for the evolution of zero-temperature 
fcrmion energy spectrum along the nodal direction in all 
cuprate superconductors. However, the thermal fluctua- 
tions are different in different materials and they destroy 
the spontaneously generated nodal gap at different tem- 
peratures. This is qualitatively in agreement with the 
ARPES data of Shen's groupi. To make quantitative 
comparison with experimental data, detailed calculations 
considering chemical structure and disorders are needed, 
which are beyond the scope of this paper. 

In order to compare our results with experiments, the 
above discussions on dynamical fermion gap generation 
should be extended to finite temperatures. The problem 
has been investigated by several authors with the results 
that above a critical temperature the chiral symmetry is 
restored and the fermions remain gapless^i. Before the 
holes are doped into the Cu-0 plane, the physical fermion 
flavor is Nf — 2. According to the results of Aitchison et 
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al. 



k B T a 
N f 



0.002, which leads to T c = 45K. From Eq. 



(19) we know that the holons shift the effective flavor to 
Nf = 3. From the result of Aitchison et al. we know 



that 



k B T c 



is about 10~ 4 , which leads to T c ~ AK. The 



nodal gap generation occurs only at temperatures below 
this critical value. When the temperature is beyond this 
critical value, thermal fluctuations destroy the nodal gap 
generation and hence the gap along the nodal direction 
closes. 



B. Competition and coexistence of 
antiferromagnetism and superconductivity 

Understanding the competition of various ground 
states of cuprate superconductors is one of the central 
problems in modern condensed matter physics. The 
cuprate superconductors in its hall-filling limit are be- 
lieved to be Mott insulators with long-range antifer- 
romagnetic order. When the doping concentration in- 
creases, long-range AF order becomes short-ranged and 
ci-wave superconductivity emerges as the ground state. 
It is interesting to build a microscopic theory to describe 
the evolution from the antiferromagnetism phase to the 
superconductivity phase with doping. 

The doping dependence of antiferromagnetism is very 
similar to that of the nodal gap observed in ARPES mea- 
surements, indicating that they might be governed by the 
same mechanism. This can be easily verified by calculat- 
ing the antiferromagnetic spin correlation function. Once 
the nodal fermion acquires a finite gap, the antiferro- 
magnetic spin correlation is greatly enhanced. Actually, 
it has been argued that the chiral symmetry breaking 
corresponds to the formation of antiferromagnetic long- 
range orde»22i2li2Si^S. The gapless spin wave excitation 
is interpreted as the massless Goldstone boson. The an- 
tiferromagnetic spin correlation is defined as 



(s + s-) = -\ 



d 3 k 



Tr[G (k)G (k+p)}, (28) 



where Go(fc) is the fermion propagator. If the fermions 
are massless, then Go(k) = ^ and we have (S + S~) = 

— yg-. At p — ► 0, (S + S~)o — > 0, and the antiferromag- 
netic correlation is heavily lost. The propagator for the 
massive fermion is 



G{k) = 



- (7 • k + «E ) 



k 2 



^0 



(29) 



where a constant mass £0 is adopted which does not 
affect the conclusion. This propagator leads to a cor- 
relation function (S + S~) that behaves like — So/27r as 
p — » and we have long-range antiferromagnetic correla- 
tion when chiral symmetry breaking takes placed. 

The gauge boson mass in the superconducting state 
is an appropriate physical quantity to describe super- 
conductivity. Thus the relationship between antiferro- 
magnetism and superconductivity can be described by 



the relationship between spontaneous nodal gap gener- 
ation and the mass of gauge boson. We have shown 
that spontaneous nodal gap generation can take place 
for doping concentration less than 8 C . This indicates 
that antiferromagnetism only exists at half-filling and 
lightly doped cuprates, in consistent with experiments. 
In most cuprate superconductors, the superfluid density 
p s is large enough to suppress antiferromagnetism once 
superconductivity emerges at the critical doping concen- 
tration, which is generally larger than 5 C . Hence, anti- 
ferromagnetism can not coexist with superconductivity 
in the bulk materials, at least in most cuprate supercon- 
ductors. However, if an external magnetic field is intro- 
duced to the superconductors and reduces the superfluid 
density p s down to below its critical value, then antifer- 
romagnetism intends to appear in the superconducting 
state. The coexistence of antiferromagnetism and super- 
conductivity is thus possible. 

Recently, intense investigation have been taken on 
the magnetic field induced local antiferromagnetic or- 
der. The external magnetic field perpendicular to the 
Cu02 plane generates Abrikosov vortices inside which 
the superfluid density is suppressed. Around the vor- 
tex cores, if the superfluid density is reduced to below 
the critical value, then spontaneous gap generation and 
hence antiferromagnetic order can be formed in the vor- 
tex stat o 4 ' 5 ' 6 ' 7 . Demler et al. have attempted to explain 
these experiments by assuming a proximity to the coexis- 
tence of spin-density wave and superconductivity^!. Such 
a coexistence can also be addressed within several other 
theoriea 32 i 33 i 34 i 35 including Zhang's SO (5) theory 3 ^ and 
the staggered flux approach proposed by Lee and Wen 3 ^. 
In this paper we use spin-charge separation and dynam- 
ical chiral symmetry breaking to account for the compe- 
tition and coexistence of antiferromagnetism and super- 
conductivity. Our approach emphasizes on the similarity 
of this competition and coexistence with the transport 
behavior and single particle spectrum properties. More 
quantitative calculation will be carried out in the future 
in order to produce the detailed experimental data. 



C. Low temperature heat transport behavior 

The low temperature transport behavior of supercon- 
ductors is controlled by the gap symmetry since it de- 
termines the type of low-energy excitations. For conven- 
tional s-wave superconductors, the low-energy density of 
states N s (uj) damps rapidly as temperature decreases, 
i.e., N s (uj) — for |o;| < Ao with Ao the quantity 
of energy gap. However, the situation is quite differ- 
ent for d-wave superconductors due to the presence gap 
nodes. In the absence of impurities, the density of states 
is Nd {oj) ~ uj. Due to the scattering of impurities a finite 
density of quasiparticle states exist, i.e., Nd (0) is finite 3 ^. 
Lee^ 7 . has investigated the quasiparticle transport of d- 
wave superconductors and found that it is independent 
of impurity intensity as a result of the competition be- 
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tween the growth of quasiparticle density of states and 
the reduction of mean free path. In contrast to s-wave su- 
perconductors where thermal conductivity vanishes due 
to the absence of mobile fcrmions at T — > 0, a linear term 
for thermal conductivity appears in d-wave superconduc- 
tors at T — > 0. Taillefer et al. for the first time observed 
such a linear term for thermal conductivity in optimally 
doped cuprate YBa2Cu306.s- 8 -. A similar linear term is 
also observed- in optimally doped Bi 2 Sr 2 CaCu208- Re- 
cently, extensive heat transport measurements have been 
performed in La2_ a: Sr a; Cu04 and YBa2Cu30j, in a wide 
range of doping concentration from overdoping to very 
underdoping. A linear term for thermal conductivity is 
observed through the whole superconducting regio n 11 ! 12 . 
Decreasing the doping concentration drives the cuprate 
to the proximity to a critical point at which sponta- 
neous gap generation occurs. As a consequence, the 
thermal conductivity should decrease down to zero at 
low doping. In the spirit of our mechanism, the gap- 
less nodal quasiparticles in the superconducting state are 
stable against gauge fluctuations because the gauge bo- 
son becomes massive and hence can not generate a finite 
gap for the gapless nodal fermions. For doping concen- 
tration below S Cl the nodal quasiparticles acquire a fi- 
nite gap which changes the node structure of d-wave gap 
symmetry and suppresses the appearance of low-energy 
fcrmions. Therefore, there is no linear term for the ther- 
mal conductivity at very low doping region. 

Recently the dependence of thermal transport on mag- 
netic field has been investigated. Naively, magnetic field 
should drive the thermal conductivity to increase since it 
breaks the Cooper pairs and hence quasiparticle density 
increases with magnetic field. This is true for many su- 
perconductors including conventional s-wave supercon- 
ductors and overdoped cuprates. However, this is not 
true for underdoped d-wave cuprate superconductors, in 
which low temperature thermal conductivity decreases 
with magnetic field. The crucial reason for this differ- 
ence is that spin-charge separation caused by strong elec- 
tron correlation takes place in underdoped cuprates but 
not in overdoped cuprates and conventional BCS super- 
conductors. When spin and charge degrees of freedom 
are separated, the holon condensation is suppressed by 
strong magnetic fields and consequently the mass of in- 
ternal gauge boson decreases down to zero with increas- 
ing magnetic field. On the contrary, the spinon pairs are 
stable against the external magnetic fields. However, the 
nodal quasiparticle is affected by the change of internal 
gauge boson mass. When superconductivity is destroyed 
by external magnetic fields at the mass gap of inter- 
nal gauge boson vanishes. As a result, the gapless nodal 
fcrmions acquire a finite gap. Below this gap the density 
of states of quasiparticles is zero, preventing the appear- 
ance of low energy fermion excitations. Thus we now 
see that, while magnetic filed generates fermion quasi- 
particles in many ordinary superconductors, the d-wave 
underdoped cuprate superconductors has a rather pecu- 
liar property that the magnetic field reduces low-energy 



fermions due to the spin-charge separation and sponta- 
neous nodal gap generation. 

When superconductivity is completely suppressed by 
magnetic field, heat is transported only by bosons includ- 
ing spin wave, holons and phonons at low temperatures. 
These bosons can only contribute a T 3 term to thermal 
conductivity at low temperature, which can be nearly 
neglected. Therefore, according to our mechanism there 
should be a thermal metal-to-insulator transition upon 
going from the superconducting state to the field-induced 
normal state of underdoped cuprates. This phenomenon 
has been observed in underdoped I^-^Sr^CuC^ re- 
cently in heat measurements^. In contrast to the ther- 
mal insulator property of the field-induced normal state, 
the holons can move freely giving rise an metal-like 
charge transport behavior. Thus we expect that the 
Wiedemann-Franz law, which gives a universal relation- 
ship between the thermal conductivity k and the elec- 
trical conductivity a as = 3- (^f ) with ks the 
Boltzmann's constant, should no longer hold in this re- 
gion. Experiments 3 ^ have provided certain evidence sup- 
porting the breakdown of this law in the electron-doped 
cuprate Pr 2 _ a; Ce 2 ;Cu04-j,. 



IV. SUMMARY 

One remarkable property of high temperature super- 
conductors is its d-wave gap symmetry. Due to the pres- 
ence of gap nodes, there is an amount of low-energy 
quasiparticles which play a crucial role in determining 
the low-temperature behavior of cuprate superconduc- 
tors. For example, the thermally excited nodal quasi- 
particles can efficiently destroy the superfluid density 28 . 
While at low temperatures, the nodal quasiparticles con- 
tribute a finite thermal conductivity that is independent 
of impurity concentrations. The d-wave gap in the single- 
particle spectrum exists not only in the superconducting 
state but also in the normal state of underdoped cuprates. 
The gapless nodal fermion excitations couples to a strong 
gauge field and can acquire a dynamically generated gap. 
When this happens, the gap nodes are removed and all 
excitations are gapped. The spontaneous nodal gap gen- 
eration modifies the picture of low-temperature physics 
to a new one, in which no free fermions can be found and 
a long-range antiferromagnetic order is formed. 

The spontaneous generation of nodal gap depends on 
the gauge boson. If the gauge boson is massless, then 
the gap generation can always take place. When super- 
conductivity emerges, the gauge boson becomes massive. 
There is a critical value for the gauge boson mass, only 
below which could spontaneous gap generation for nodal 
fermions take place. The high temperature superconduc- 
tors have a peculiarity that the gauge boson mass is pro- 
portional to the doping concentration. Thus spontaneous 
nodal gap generation should exist at low doping concen- 
trations while vanishes as doping increases. However, 
when superconductivity is suppressed by some external 
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means such as external magnetic field, spontaneous nodal 
gap generation occurs and correspondingly low temper- 
ature nodal quasiparticles disappear. Hence there is a 
thermal metal-insulator transition from the supercon- 
ducting ground state to the field-induced normal ground 
state upon increasing magnetic field. On the other hand, 
spontaneous nodal gap generation corresponds to the for- 
mation of antiferromagnetic long-range order. The fate 
of antiferromagnetism hence can be described by the re- 
lationship of spontaneous nodal gap generation with su- 
perfluid density or doping concentration. 



ACKNOWLEDGEMENT 



G.-Z. Liu would like to thank Prof. X.-H. Chen for 
helpful discussions. G.-Z. Liu is supported by the Re- 
inforcing Education Project Towards 21th Century. G. 
Cheng is supported by National Science Foundation in 
China No.10175058 and No.10231050. 



1 K. M. Shen, T. Yoshida, D. H. Lu, F. Ronning, N. P. 
Armitage, W. S. Lee, X. J. Zhou, A. Damascelli, D. L. 
Feng, N. J. C. Ingle, H. Eisaki, Y. Kohsaka, H. Tak- 
agi, T. Kakeshita, S. Uchida, P. K. Mang, M. Greven, 
Y. Onose, Y. Taguchi, Y. Tokura, S. Komiya, Y. Ando, 
M. Azuma, M. Takano, A. Fujimori, and Z.-X. Shen, 
cond-mat/0312270 

2 C. C. Tsuei and J. R. Kirtley, Rev. Mod. Phys. 72, 969 

(2000) . 

A. Damascelli, Z. Hussain, and Z.-X. Shen, Rev. Mod. 
Phys. 75, 473 (2003). 

4 B. Lake, G. Aeppli, K. N. Clausen, D. F. McMorrow, K. 
Lefmann, N. E. Hussey, N. Mangkorntong, M. Nohara, H. 
Takagi, T. E. Mason, and A. Schroder, Science 291, 1759 

(2001) ; B. Lake, H. M. Ronnow, N. B. Christensen, G. 
Aeppli, K. Lefmann, D. F. McMorrow, P. Vorderwisch, 
P. Smeibidl, N. Mangkarntong, T. Sasagawa, M. Nohara, 
H. Takagi, and T. E. Mason, Nature (London) 415, 299 

(2002) . 

5 V. F. Mitrovic, E. E. Sigmund, M. Eschrig, H. N. Bach- 
man, W. P. Halperin, A. P. Reyes, P. Kuhns, and W. G 
Moulton, Nature (London) 413, 501 (2001). 

6 J. E. Hoffman, E. W. Hudson, K. M. Lang, V. Madhavan, 
H. Eisaki, S. Uchida, and J. C. Davis, Science 295, 466 
(2002). 

7 B. Khaykovich, Y. S. Lee, R. W. Erwin, S.-H. Lee, S. Waki- 
moto, K. J. Thomas, M. A. Kastner, and R. J. Birgeneau, 
Phys. Rev. B 66, 014528 (2002). 

8 L. Taillefer, B. Lussier, R. Gagnon, K. Behnia, and H. 
Aubin, Phys. Rev. Lett. 79, 483 (1997). 

9 M. Chiao, R. W. Hill, Ch. Lupien, L. Taillefer, P. Lambert, 
R. Gagnon, and P. Fournier, Phys. Rev. B 62, 3554 (2000). 

10 D. G. Hawthorn, R. W. Hill, C. Proust, F. Ronning, M. 
Sutherland, E. Boaknin, C. Lupien, M. A. Tanatar, J. 
Paglione, S. Wakimoto, H. Zhang, L. Taillefer, T. Kimura, 
M. Nohara, H. Takagi, and N. E. Hussey, Phys. Rev. Lett. 
90, 197004 (2003). 

11 J. Takeya, Y. Ando, S. Komiya, and X. F. Sun, Phys. Rev. 
Lett. 88, 077001 (2003). 

12 M. Sutherland, D. G. Hawthorn, R. W. Hill, F. Ronning, S. 
Wakimoto, H. Zhang, C. Proust, E. Boaknin, C. Lupien, L. 
Taillefer, R. Liang, D. A. Bonn, W. N. Hardy, R. Gagnon, 
N. E. Hussey, T. Kimura, M. Nohara, and H. Takagi, Phys. 
Rev. B 67, 174520 (2003). 

13 P. W. Anderson, Science 235, 1196 (1987). 

14 F. C. Zhang and T. M. Rice, Phys. Rev. B 37, 3759 (1988). 

15 G. Baskaran, Z. Zou, and P. W. Anderson, Solid State 



Commun. 63, 973 (1987). 

16 G. Kotliar and J. Liu, Phys. Rev. B 38, 5142 (1988). 

17 I. Affleck and J. B. Marston, Phys. Rev. B 37, 3774 (1988); 
J. B. Marston and I. Affleck, Phys. Rev. B 38, 11538 
(1989). 

18 L. B. Ioffe and A. I. Larkin, Phys. Rev. B 39, 8988 (1989). 

19 X.-G. Wen and P. A. Lee, Phys. Rev. Lett. 76, 503 (1996); 
P. A. Lee, N. Nagaosa, T. K. Ng, and X.-G. Wen, Phys. 
Rev. B 57, 6003 (1998). 

20 D. H. Kim and P. A. Lee, Ann. Phys. (N.Y.) 272, 130 
(1999). 

21 G.-Z. Liu and G. Cheng, Phys. Lett. B 510, 320 (2001). 

22 T. Appelquist, D. Nash, and L. C. R. Wijewardhana, Phys. 
Rev. Lett. 60, 2575 (1988). 

23 Some recent studies on QED3 have arrived results against 
the conclusion that the critival fermion flavor is about 
3.3. Numerical simulations on lattice QED3 performed by 
Hands et al. found no decisive signal for chiral symmetry 
breaking for N > 2^. It was suggested by Gusynin and 
Reendera— that this result is owing to the introduction of 
large infrared cutoff in lattice studies. Several years ago, 
Appelquist et a/4- argued that chiral symmetry of QED3 
can not be broken if the flavor is larger than 3/2 based on a 
conjectured constraint on the infrared structure of asymp- 
totically free gauge theories. However, this result is doubful 
because the reliability of this constraint is not verified, as 
mentioned by Appelquist et al. in a very recent review—. 
On the other hand, the original result in Ref. (22) has 
been confirmed by extensive investigations—*— and now is 
widely accepted. 

24 G.-Z. Liu and G. Cheng, Phys. Rev. B 66, 100505 (2002). 

25 G.-Z. Liu and G. Cheng, Phys. Rev. B 65, 132513 (2002). 

26 G.-Z. Liu and G. Cheng, Phys. Rev. D 67, 065010 (2003); 
G.-Z. Liu and G. Cheng, cond-mat/0208061 

27 J. Orenstein, G. A. Thomas, A. J. Millis, S. L. Cooper, 

D. H. Rapkine, T. Timusk, L. F. Schneemeyer, and J. V. 
Waszczak, Phys. Rev. B 42, 6342 (1990). 

28 P. A. Lee and X.-G. Wen, Phys. Rev. Lett. 78, 4111 (1997). 

29 I. J. R. Aitchison, N. Dorey, M. Klein-Kreisler, and N. 

E. Mavromatos, Phys. Lett. B 294, 91 (1992); I. J. R. 
Aitchison and M. Klein-Kreisler, Phys. Rev. D 50, 1068 
(1994). 

30 Z. Tesanovic, O. Vafek, and M. Franz, Phys. Rev. B 65, 
180511 (2002); I. F. Herbut, Phys. Rev. Lett. 88, 047006 
(2002). 

31 E. Dernier, S. Sachdev, and Y. Zhang, Phys. Rev. Lett. 87, 
067202 (2001). 



10 



S.-C. Zhang, Science 275, 1089 (1997); D. P. Arovas, A. 
J. Berlinsky, C. Kallin, and S.-C. Zhang, Phys. Rev. Lett. 
79, 2871 (1997). 

J.-I. Kishine, P. A. Lee, and X.-G. Wen, Phys. Rev. Lett. 
86, 5365 (2001). 

M. Franz and Z. Tesanovic, Phys. Rev. Lett. 88, 257005 
(2002). 

J.-X. Zhu and C. S. Ting, Phys. Rev. Lett. 87, 147002 
(2001). 

A. C. Durst and P. A. Lee, Phys. Rev. B 62, 1270 (2000). 
P. A. Lee, Phys. Rev. Lett. 71, 1887 (1993). 
R. W. Hill, C. Proust, L. Taillefer, P. Fournier, and R. L. 
Greene, Nature (London) 414, 711 (2001). 



S. Hands, J. Kogut, and K. Strouthos, Nucl. Phys. B 645, 
321 (2002). 

V. P. Gusynin and M. Reenders, Phys. Rev. D 68, 025017 
(2003). 

T. Appelquist, A. Cohen, and M. Schmaltz, Phys. Rev. D 
60, 045003 (1999). 

T. Appelquist and L. C. R. Wijewardhana, 
hep-ph/0403250 

D. Nash, Phys. Rev. Lett. 62. 3024 (1989). 

P. Maris, Phys. Rev. D 54, 4049 (1996); V. Gusynin, A. 

Hams, and M. Reenders, Phys. Rev. D 63, 045025 (2001). 



